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Abstract 

Superstatistics and Tsallis statistics in statistical mechanics is given an interpretation in terms of 
Bayesian statistical analysis. Subsequently superstatistics is extended by replacing each component of the 
conditional and marginal densities by Mathai's pathway model and further both components are replaced by 
Mathai's pathway models. This produces a wide class of mathematically and statistically interesting func- 
tions for prospective applications in statistical physics. It is pointed out that the hnal integral is a particular 
case of a general class of integrals introduced by the authors earlier. Those integrals are also connected to 
Kraetzel integrals in applied analysis, inverse Gaussian densities in stochastic processes, reaction rate inte- 
grals in the theory of nuclear astrophysics and Tsallis statistics in nonextensive statistical mechanics. The 
final results are obtained in terms of Fox's H-function. Matrix variate analogue of one significant specific 
case is also pointed out. 
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1. Introduction 

Nonequilibrium systems in various areas of science, see Ebeling and Sokolov (2005) and Uffink (2007), 
exhibit spatio-temporal dynamics that is inhomogeneous and can be described by a superposition of several 
statistics on different scales, in short called a superstatistics approach in the physical literature, see Beck 
(2006). Essential for this superstatistics approach is the existence of sufficient time scale separation between 
two relevant dynamics within the nonequilibrium system. There must be some intensive parameter that 
fluctuates on a larger time scale than the typical relaxation time of the local dynamics. In a statistical 
thermodynamic setting such a parameter can be interpreted as a local inverse temperature of the system, 
but much broader interpretations are possible dependend on the nonequilibrium system under consideration. 
The stationary distributions of superstatistical systems, obtained by averaging over all relevant fluctuating 
parameters exhibit non-Gaussian behaviour with fat tails, which can be a power law, or a streched exponen- 
tial, or other functional forms like Mittag-Leffler function or Fox H-function. The relevant superstatistical 
paremeter can be an effective parameter in a stochastic differential equation or a local variance parameter 
extracted from a time series. Currently, in the physical literature, superstatistics [Beck (2005)] and Tsal- 
lis statistics [Tsallis (1988)] are the two preferred models to describe nonequilibrium systems in a sense of 
generalyzing the well-established Boltzmann-Gibbs statistical mechanics. 

In this paper we shall develop a Bayesian approach to superstatistics and Tsallis statistics in terms of 
the recently introduced pathway model of Mathai (2005). It will be shown that Mathai's pathway model 
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comprises the distributions of superstatistical systems as well as Tsallis statistics in a unified manner. 



2. Bayesian Statistical Analysis 

Let us start with the standard Bayesian analysis problem. Consider a positive real scalar random 
variable x and a parameter 9. Let the conditional density of x at a given value of 9 be denoted by f(x\8) 
and the marginal density of 9 by g(9) respectively. For simplicity, let us take both f{x\9) and g(9) in the 
generalized gamma family of functions. Let 



for 9 > 0, 5 > 0, x > 0, 7 > and f{x\9) = elsewhere. Let 
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for b > 0, 9 > 0, S > 0, p > and 5(6*) = elsewhere. Then the unconditional density of x is given by the 
following: 
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for x > 0, 7 > 0, p > 0, S > 0, 6 > 0. 

Theorem 1. Le£ i/ie conditional density of x given 9, that is, f(x\9), and the marginal density of 9 be as 
in (1) and (2) respectively. Then the unconditional density of x, denoted by f x (x), is given by (3). 

In a physical system the parameter 9 in (1) may represent temperature so that the density f{x\9) may 
represent the production of the item x at a fixed temperature 9 or at a given value of 9. Then the marginal 
density of 9 in (2) may represent the temperature distribution. What is the distribution of the production 
of x over all temperature variations or averaged over the density of 91 This is the unconditional density of x 
given in (3) for the specific densities in (1) and (2). Since a density such as the one in (2) is superimposed over 
the density such as the one in (1), the resulting density in (3) is called superstatistics. Various interpretations 
of x and 9 in different physical systems may be seen in the original paper on superstatistics is that of Beck 
and Cohen (2003). Equation (3) for 5 = 1, 7 = 1, b = q > 1 and = is Tsallis statistics of 

non-extensive statistical mechanics, see Tsallis (1988). The Bayes' density of 9 is 

= (4, 

Jx\X) 

which is available from (1),(2) and (3). The Bayes' estimate of 9 is the conditional expectation of 9 at given 
value of x, denoted by E{9\x), where E denotes statistical expectation. This is available from the conditional 
density of 9, given x, in (4). Thus, Bayes' procedure, superstatistics and Tsallis statistics are all connected 
together as shown in (1) to (4). 

3. Extension of Bayes' Procedure 

We can extend the discussion in (1) to (4) by using the pathway model of Mathai (2005). For example, 
let us replace g(9) in (2) by a pathway density, namely, 
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for /3 > 1, 5 > 0, r?i > 0, 6 > 0, p > 0, 9 > 0, ^ - § > and Pi(0) = elsewhere. When -> 1+ in 
(5) the model in (5) goes to <?(0) in (2) with 6 replaced by br\\. The model in (5) consists of three different 
functional forms. For /? < 1 we may write 



[1 + b(0 - l)0 d p~ = [f - 6(1 - I3)e 6 ] — . 

The right side remains positive in the finite range f — 6(1 — /3)9 S > or < < [6(1 — f3)]~i . Then for j3 < 1 
the density in (5) changes to the form 



^(i -^r^ + i + f) 
r(f)r(^ + i) 



W) = FTZwj^n -^"Ml - 6(1 - 0)^- (6) 



Note that for /3 > 1, the density in (5) stays in the generalized type-2 beta family of densities and for < 1 
the density in (6) belongs to the generalized type-1 beta family of densities. When f3 — > 1, either from the 
left or from the right, P 2 {9) and P\{9) will go to Pz{9), where 

Ps(0) = S -^9^e-^\ (7) 

for 6 > 0, ryi > 0, S > 0, p > 0, > 0, which is the generalized gamma family of densities. Thus, the 
parameter [3 creates a path of going from the generalized type-1 beta family to generalized type-2 beta 
family to generalized gamma family of functions. It is shown in Mathai and Haubold (2007) that almost all 
statistical densities in current use in different disciplines are available from (5). In a physical system, if (7) 
is the stable situation then the unstable neighborhoods are covered by (5) and (6) and the paths in between, 
the path described by the parameter /3. Note that the model in (6) is a model with the right tail cut off. 
The movement of /3 will provide thicker or thinner-tailed distributions which will also be helpful in a model 
building situation where one may be looking for a thicker or thinner-tailed distribution as an appropriate fit. 

It is not difficult to show that the normalizing constants in (5) and (6) reduce to the normalizing constant 
in (7). This can be seen by expanding the gammas with the help of Stirling's formula, which is given by 

T(z + a) w V2^z z+a -h~ z 

for \z\ -4- oo and a a bounded quantity. When /? — > 1_ we have — >■ oo and when f3 -4- 1 + , -4- oo. 
Hence take z = or z — -^y as the case may be, expand and simplify to see that the normalizing 
constants in (5) and (6) go to the constant in (7). 

Instead of g(0) of (2) if we use the extended density Pi(0) of (5) then the unconditional density f x {x) 
in (3) has the following form: 
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The integral part in (8) can be evaluated by using Mellin convolution property by observing that the integrand 
is of the form 9hx(x9)h 2 {9) for some functions h\ and hi. Let 



h x {9) = ciflT-V " and h 2 {9) = c 2 9 p - 1 [l + b(J3 - 1)0' 5 ]" 



3 



so that hi and h 2 can create statistical densities for appropriate values of the normalizing constants c\ and 
c 2 . Therefore 
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Hence the Mellin transform of the left side is the product of the Mcllin transforms of the right side with 
parameters s for hi and 2 — s for h 2 - Writing in terms of statistical expectations, 
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for — s + 1) > 0, 5£(^rj — p ~^ +1 ) > 0. Therefore the inverse Mcllin transform is the following: 
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where H is the H-function, see for example, Mathai, Saxena, and Haubold (2010) and Srivastava, Gupta, 
and Goyal (1982). 

Theorem 2. When the conditional density of x given 9 is of the form f(x\9) in (1), when the marginal 
density of 9 is of the form of Pi (6) in (5) then the unconditional density of x, denoted by f x (x) is given by 
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/or 6 > 0, /? > 1, 7 > 1, p > 0, 5 > 0, ^ - £+1 > 0. 

4. Pathway Extension of the Conditional Density 

Let us consider an extended form of the conditional density f{x\9) and the original form of the marginal 
density as in (2). Let the extended conditional density be of the pathway model type-2, given by, 



S[a(a- l)]^r(^L) . . „ 



(10) 



for a > 1, 5 > 0, r\ > 0, 7 > 0, > 0, x > and f(x\0) — elsewhere. Let the marginal density of be as 
in (2). Then the unconditional density of x, again denoted by f x {x), is given by the following: 



f x (x) = / f(x\0)g(0)d0 
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The integral part can be written as follows: 
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where 

hi(xi) = c\xl~ l [l + a(a - l)^] - ^ 

and 

h 2 {x 2 ) = c 2 x^ 1 e- bx K 

Then c\ and C2 can act as normalizing constants so that hi and h 2 are statistical densities. Then, writing 
the Mellin transform in terms of expected values, we have the following: 
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Hence the inverse Mellin transform of E^ -1 )^^" 8 ) is given by 
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Theorem 3. Lei i/ie conditional density of x at given be given by (10) and the marginal density of be 
given by (2). Then the unconditional density of x, denoted by f x (x), is given by 
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5. General Pathway Model for Unconditional Densities 

Now we consider a more general case where f(x\0) is given in (10) and g(0) is given in (5). Then what 
will be the unconditional density of xl From the procedure so far, it is clear that the unconditional density 
of x, denoted by f x (x),is given by the following: 
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for > 1, a > 1, 5 > 0, t? > 0, 771 > 0, b > 0, a > 0, p > 0, 6 > 0, or > 0, ^_ - J > 0, ^ - § > 0. 
Consider the integral part. 
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where 
and 



h 1 (x 1 )=c 1 xj 1 [l + a(a- l)^]"^ 
h 2 (x 2 ) = c 2 x p 2 - 1 [l + b(0 - l)x s 2 ] -7»=r. 



When ci and c 2 are normalizing constants then h\(x\) and ^2(^2) can act as statistical densities. 
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Then the integral is available from the inverse Mellin transform which can again be written in terms of a 
H-function. 



Theorem 4. Let the conditional density of x given 8 be given in (10) and the marginal density of be 
given in (5). Then the unconditional density, denoted by f x (x), is given by 
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Remark 1. The integrand in this general case, giving rise to Theorem 4, is in the form of a product of two 
pathway models. This, in fact is a special case of a versatile integral considered in Mathai (2007). This will 
be briefly discussed in the next section. 



6. A Versatile Integral 

One form of the versatile integral is the following: 

POO 

f(z2\zi)= / x^ll + zfia-^x^-^il + z^-^x-T^dx (15) 
Jo 

where the parameters a and [3 can be independently less than one, greater than one and going to one. Also 5 
and p can be independently positive or negative, all are assumed to be real quantities for convenience. Thus 
the integral in (15) covers a large spectrum of integrals. The integrand in (15) is nothing but a product of 
two pathway models in the real scalar case. When a — > 1 and (3 — > 1 we have several important integrals in 
different fields: (1): Kratzel integral in applied analysis is a particular case in this situation with 5 = 1 and 
p = 1. (2): When p = \ and <5 = 1 we have the reaction rate probability integral in the theory of nuclear 
astrophysics discussed in a series of papers, see for example, Mathai and Haubold (1988). (3): When p = 1 
one has the inverse Gaussian density, which appears very frequently in time series, stochastic processes and 
statistical distribution theory. Generalizations of all these practical situations are given by (15) for various 
values of a, (3, S, p. 

7. Matrix-variate Analogues 

All cases of (1) to (15) cannot be given matrix-variate analogues due to restrictions on Jacobians of 
matrix transformations, see for example Mathai (1997). One situation is straightforward. Consider real 
positive definite p x p matrices X and Y having densities belonging to the matrix-variate gamma families. 
Let the conditional density of Y, given X, be given by 

\x\ — 

f(Y\X) = l —^—\XY\ a - E i 1 e- t ^ XY) , X>0, Y > 
r p (a) 

and 

g{ X) = ^\Xr^e-^\ B>0 

where |(-)| denotes the determinant of (•), tr(-) the trace of (•) and the real matrix-variate gamma function 
is given as 

T p (a) = n^T(a)T(a l)..T(a - ^), R(a) > P -^. 
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The standard notation A > is used for the matrix A being real symmetric, and further, positive definite. 
Then the unconditional density of Y, denoted by fy(Y), is given by the following: 

MY) = [ f(Y\X)g(X)dX 
Jx 

= - } B X, , 1 lYr-^ixr+e-^e-^+^dx 

T p (a)T p (P) J x 

-iv&^* |fl+11 ~ c " Br ' (0+ffl 

= WTOW^i' + B -^-" +>> - r>°.*>°- ('«) 

The transformation used in the integration is U = (B + Y)^X(B + Y)^ dX = \B + Y\~ E i~dU where, 
for example, dU denotes the wedge product of the p(p + l)/2 differentials in the elements of U and (•) 2 
denotes the positive definite square root of the positive definite matrix (•). From (16) it is clear that the 
unconditional density of Y is a matrix- variate type-2 beta density. This is the matrix- variate analogue of 
the superstatistics in the case of conditional and marginal densities belonging to the gamma type densities. 
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